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Abstract: When using dimensional regularization/reduction the e-dimensional numerator 
of the 1-loop Feynman diagrams gives rise to rational contributions. I list the set of 
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1. Introduction 



|, ^, 10, 11, 12, for computing 1-loop corrections 

that, directly or indirectly, has allowed an impressing 



New techniques |, |, |, |, |, § g 
led to a NLO revolution ^ |l§, 
improvement in our ability to predict physical observables at the NLO accuracy [|17| , |18 
ll, p|, M M |2|, M |5l g 1^, ^, p9|, |30, |3l| 



a PI, 35, 3e, 



. Basically all 

new methods need a special treatment of the contributions that are not proportional to 
the scalar 1-loop functions, the so called rational terms [ 



39, 40, 41 1 . That is achieved, 



in Unitarity and Generalized Unitarity methods, by computing the entire amplitude in 
different numbers of space-time dimensions |42], or via bootstrapping techniques p3|, E3|, 



or through d-dimensional cuts [45|. The Ossola-Papadopoulos-Pittau (OPP) approach of 
reference requires, instead, the computation (once for all for the theory at hand) of a 
special set of tree level Feynman rules ||4^, |4^, up to 4-point interactions ^ . 

When using dimensional regularization/reduction, the origin of such terms lies in the 
e-dimensional numerator of the 1-loop Feynman diagrams ^. To be more specific, let us 
consider the general expression for the integrand of a generic m-point one-loop (sub-) 
amplitude 



N{q) 



DqDi 



D. 



m— 1 



Di = {q+ Pi 



[i.D 



^The contributions to higher-point functions vanish because of UV finiteness. 

^Another contribution, which is however directly linkable to the cut-constructible part of the amplitude 
5l|l, is generated by the e-dimensional 1-loop denominators. 
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where q is the integration momentum and 



-2 2,-2 

1 = q +q 



2 2 

7 -/^ 



(1.2) 



In the previous expressions and in all the following ones, a bar denotes objects living in 
d = 4 + e dimensions, whereas a tilde represents e-dimensional quantities. 

The numerator function N(q) can be split into a 4-dimensional plus an e-dimensional 

part 



Niq) = Niq) + Nif,^q,e). 



(1.3) 



N{q) lives in 4 dimensions, while N{fj? ,q,e), which originates from the splitting of d- 
dimensional objects 



Q = q + q, 

gives rise to a rational piece called R2 in the OPP language: 



HV (27r)4 



with 



DqDi ■ ■ ■ An-l 



d'^q / d'^jL . 



(1.4) 



(1.5) 



(1.6) 



R2 has a pure ultraviolet origin ^] , so that eq. (^]^) also holds when infrared /collinear 
divergences are present in the loop integrals. It can be shown [39, 42| that iV(/u^,g, e) is 



polynomial in /u^ and at most linear in e, and that the e dependence can be reabsorbed 
in the regularization scheme ||52|. Therefore, beside eq. (|1.5| ), which defines R2 in the 't 
Hooft-Veltman (HV) scheme, a Four Dimensional Helicity scheme (FDH) can be used in 
which the e dependence in the numerator function is discarded before integration 



1 



FDH {2ttY 



0) 



DqDi ■ ■ ■ Z)^_ 



m— 1 



(1.7) 



As for the virtual part of the NLO corrections, FDH is equivalent to Dimensional Reduc- 
tion [53|. 

It is clear that explicitly using the rules in eq. ( |1.4| ) allows an analytic extraction, 
Feynman diagram by Feynman diagram, of the coefficients of the various powers of fj? 
and e. For example, the GoSam [^, 55 1 approach achieves this on the fly, by linking 
to algebraical manipulation programs providing the necessary algebra when building the 
amplitude. The case of QCD is particularly simple, and computations based on a standard 



Passarino-Veltman |56| decomposition are relatively easy [g^, 58]. In addition, for gluonic 



amplitudes, a super-symmetric decomposition relates the contribution of the rational terms 
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to a scalar massive gluon running in the loop which can also be computed by using 
massive Cachazo-Svrcek-Witten Feynman rules |6^, |62| . 

All the approaches mentioned so far have advantages and drawbacks. It is however 
beyond doubt that it would be desirable to have a way to compute R2 independent of the 
theory at hand, four dimensional and not requiring the use of analytical manipulations. 
In this paper, I provide such a method in the form of primary Feynman rules, which 
reproduce the polynomial dependence on of the integrand in eq. ( |1.7| ) ^ when working in 
the renormalizable gauge. As will be explained in the next section, such rules are uniquely 
determined by reading the original propagators and vertices of the theory, and solely depend 
on their Lorentz structure. They can therefore be used as any other Feynman rules by 
programs such as MadLoop [12|, HELAC-NLO [^], FeynRules [63] or within a methods such 



as Open Loops [^] to automatically generate R2. 

This paper is organized as follows: in section ^ I give an introductory discussion and fix 
my conventions. Section |^ introduces vectors, scalars and their mutual interactions. The 
fermionic case is more delicate, and it is discussed in section ^. In section |5| the generic 
formulae are specialized to the QCD case and, in section ^, I draw my conclusions. Further 
details are reported in appendix 0. 

2. Preliminaries and conventions 

I am interested in reconstructing the polynomial dependence on /i^ of the e-dimensional 
numerator in eq. ([1.7]). In a renormalizable gauge the relevant expansion is 

2 

N{f^',q,e = 0) = Y,il^ycj{q). (2.1) 

i=i 

As already stated in section |l|, an explicit dependence on e, such as that implied by eq. (|L^), 
can always be reproduced by a change of regularization scheme. The translation rules 
relevant for QCD and for the electroweak standard model are given in appendix |A[ 
By inserting eq. ( 2A ) into eq. (^^) the relevant integrals are of the kind 



^d-i^y^i::^^ ^.^j^ m<4,0<j<2,r<2, (2.2) 
(27r)4y ^DoDi---Dm-i " ' - ' - ' ^ ' 

which give a non- vanishing contribution, in the limit e — t- 0, only when 4 + 2j + r — 2m > 0. 
They are computed in |5l[ ]. 

Powers of jj,"^ in eq. ( |2.lD are generated by the contraction of the integration momen- 
tum with itself only in the presence of vectors and fermions, which are indeed the only 
particles bringing the integration momentum in the numerator. For each vector field V^, 
I therefore introduce a scalar field V, whose propagator corresponds to the propagation of 
its e-dimensional components. In the same way, for any fermion F, I introduce a fermionic 
field F and its corresponding propagator. In the following, V and F are called e-particles 



^The conversion to the HV scheme of eq. (1.5) is presented in appendix ^ 
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and their propagators are graphically represented by a dashed line and a dashed arrow, 
respectively. 

e-particles interact among themselves and with the particles of the original theory, 
bringing an explicit dependence on fi and are only allowed to circulate in the loop. I call 
e-vertices the special vertices involving e-particles, and e-diagrams the Feynman diagrams 
that contain e-vertices. Given the fact that the only possible mechanisms for generating 
powers of /i^ are those illustrated in eq. ( |1.4| ) , it is easy to determine the e-vertices by simply 
looking at the Feynman rules of the original theory. In particular, the Lorentz and color 
structures are completely dictated by the nature of the particles entering the e-vertices, so 
that, in general, only the relative phases of /i and/or have to be determined. I fixed 
them by explicitly writing down all possible classes of 2,3, and 4-point Feynman diagrams 
involving scalars, vectors and fermions and by requiring that the e-vertices reproduce the 
results obtained with an explicit calculation of R2. It should be possible to determine the 
e-vertices directly from the original Lagrangian, by splitting it into 4 and e dimensional 
parts, but I did not choose such an approach. In the following two sections, I list all 
possibilities involving vectors, scalars and fermions. I do it in a completely generic way, in 
the sense that any particular model can be implemented just by giving a specific value to 
the constants in front of the listed Lorentz structures. For example, the electroweak model 
is obtained by fixing them according to reference ||6^ and the minimal supersymmetric 
standard model by using the rules in []6^ ]. 



3. Vectors, scalars and their interactions 

The propagator of a scalar e-particle associated with a vector field is given in fig. ||. From 
the original three-vector vertex, two corresponding e-vertices are derived, as illustrated 
in fig. ^. To determine the sign of the first e- vertex, one should keep track of the fiow of the 
loop momentum q. The original four-vector vertex gives rise to the two e-vertices of fig. |3|, 
while two-vector-one-scalar, one-vector-two-scalar and two-vector-two-scalar vertices gen- 
erate just one e- vertex each, as shown in figs. ^ to ^. 

Although the e-vertices can be used in any stage of the calculation to compute e- 
diagrams reproducing the coefficients Cj{q) of eq. ( |2.1| ), it may be useful and illuminating 
to consider them in strict connection with the original Feynman diagrams. With this 
kind of reasoning, one may associate to any 1-loop Feynman diagram contributing to the 
amplitude under study, a set of e-diagrams which fully reconstruct its /i^ dependence. In 
particular, when using the same rooting for the loop momentum in each of them, this 
reconstruction even holds at the integrand level. Furthermore, given the fact that the e- 
vertices bring known powers of fi into the e-diagrams, one can easily disentangle sub-sets 
generating specific powers of /i^ in eq. {2A). A concrete example is given in fig. for two 



particular diagrams contributing to a VV — >• VV scattering, where the sum of the two 
e-diagrams in the second line yields the term C2{q) corresponding to the original box 
diagram. 
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V 



V 



Figure 1: Propagator of a vector particle (box on tlie top) and propagator of its corresponding 
scalar e-particle. 
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Figure 2: Three-vector vertex (box on the top) and its corresponding e-vertices. q represents the 
flow of the loop momentum. 



4. Fermions and their interactions with vectors and scalars 



The study of the interactions involving fermions is complicated by the presence of 75 . It is 
convenient to start from the chiral fermions of the original theory and to split the fermion 
propagator in chirality flipping and chirality preserving parts, as illustrated in the top part 
of fig. ^. The e-propagator corresponding to the e-particle associated with a fermion is 
chirality flipping, being reminiscent of the presence of a ^ in the numerator, as shown in 
the bottom part of fig. |8[ Fermions can interact with vectors and scalars with the standard 
vector-fermion-fermion and scalar-fermion-fermion vertices shown in the top parts of figs. 



and |10|, while the corresponding e-vertices are drawn in the bottom parts. Notice also that 
all vertices in fig. |^ are chirality flipping, due to the presence of 7^ in the original interaction, 
while those in fig. |l^ are chirality preserving, because of their scalar nature. Particularly 
interesting is the vertex in fig. ^ (c), which represents a VFF interaction. Although no 7 
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Figure 3: Four-vector vertex (box on the top) and its corresponding e-vertices. 
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Figure 4: Two-vector-one-scalar vertex (box on the top) and its corresponding e-vertex. 



matrix is present in that vertex, it should be considered as a chirality flipping one, because 
of its vectorial origin. 



e-diagrams are built by using the rules of figs. |8| to |10| and reading, as usual, the 
fermionic line backward starting from the arrow. After the last vertex is encountered, 
a chirality projector = i(l ± 75) should be inserted, according to the chirality of 
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Figure 5: One- vector- two-scalar vertex (box on the top) and its corresponding e- vertex, q 
represents the flow of the loop momentum. 
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Figure 6: Two-vector-two-scalar vertex (box on the top) and its corresponding e-vertex. 



the fermion entering into it. For example, the diagram in fig. 11 generates the fermionic 
structure 

^^(l)7^7a7^^^«(2) . (4.1) 
For fermionic loops, the starting point is arbitrary and should be kept fixed when summing 



over diagrams and families in order to preserve the right cancellations |47|. 

When scalars are present, the rules presented in figs. |8| to 10 define the diagrams up to 
a sign, which has to be included by hand. This is due to the anticommutation properties 
of I with the four dimensional 7 matrices. The rule for fixing the sign is as follows. For 
each pair of e-fermion lines present in a given diagram the result should be multiplied 
by (— )("3+"p)j where is the number of scalar vertices and Up the number of chirality 
preserving propagators of the kind i- 



^ between them. For example, a minus sign should 



-7- 




Figure 7: Examples of e-diagrams. In the two boxes I draw the original diagrams. The e-diagrams 
below each box reconstruct the complete iJ? dependence of the integrand of the corresponding 
diagram, provided the same rooting for the loop momentum is chosen in each of them. 



be assigned to the diagrams in fig. (a), (b) and (d) whereas a plus sign should be given 
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Figure 8: In the box on the top a fermion propagator is split in chirality flipping and chirality 
preserving parts, h = ± denotes right-handed or left-handed components and the dashed line on 
the bottom represents the propagator of the e-particle associated with a fermion. 



to that one in fig. (c) . 

A last subtlety concerns again the vertex in fig. |9| (c). The scalar e-dimensional degrees 
of freedom brought by the field V can only originate from one of the vertices introduced 
in section ^ In fact, V fields are only needed for non abelian theories because of the 
momentum dependent three- vector vertex of fig. ^ ^. Therefore, diagrams where both ends 
of a y e-particle connect to a fermion line should be discarded, such as that one given 



in fig. 13 



Finally, in the case of Major ana fermions, such as neutralinos and gluinos in SUSY 
theories, the relative sign of interfering Feynman diagrams can be determined as described 
in 1^,0. 

5. QCD 

In the case of QCD no 75 is present and splitting fermions into right-handed and left- 
handed components is no longer necessary. No additional sign needs to be inserted, since 
no scalar particles are involved in the e-fermionic vertices, although diagrams in which a 
scalar e-gluon connects 2 fermionic lines should be discarded, as explained in the previous 



section. I list the relevant QCD e-propagators and e-vertices in fig. As can be seen, 
they have exactly the same general Lorentz structure described in sections |3| and while 
terms of the original color structures split among the various contributions in a well defined 
way. 

6. Conclusions 

I presented the set of special Feynman rules allowing the reconstruction of the e-dimensional 
part of 1-loop amplitudes in theories with vectors, scalars and fermions. The rules are quite 



^In pure QED, one just needs to introduce one F field for each fermion family. 
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Figure 9: Vector-fermion-fermion vertex (box on the top) and its corresponding e-vertices. h — ± 
denotes right-handed or left-handed fermion components, ujh is a chirahty projector and q represents 
the flow of the loop momentum. 
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Figure 10: Scalar- fermion- fermion vertex (box on the top) and its corresponding e-vertices. h = ± 
denotes right-handed or left-handed fermion components and ujh is a chirality projector. 



simple, when assuming a renormalizable gauge, and easily derivable from the vertices of 
the original theory. They can be used to extract the /i^ dependence from the integrand 
of any contributing 1-loop Feynman diagram, namely the e-dimensional part generated by 
self contractions of the loop momentum. 
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Figure 11: e-diagram generating the fermionic structure 11(1)7/3707'^ a;'*u(2) 
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Figure 12: Fermionic e-diagrams. A minus sign should be assigned to the diagrams (a), (b) and 
(d), while no additional sign is required for (c). 



Figure 13: Example of fermionic e-diagram that should be discarded because both ends of the 
e-gluon connect to a fermionic line. 



The complete electroweak model can be studies by simply fixing the constants appear- 
ing in the vertices of figs. |l] to ^ and figs. |8| to IC to their standard model values, while the 



interactions relevant for QCD are explicitly listed in fig. |14. 

A four dimensional helicity scheme is used in this work, but simple translation rules 
to the 't Hooft Veltman scheme are collected in an appendix. 

SUSY and BSM theories sharing the same Lorentz structures studied in this paper can 
be treated in the same way. 

The special vertices presented here may also be considered as a practical tool to de- 
termine the counter-terms needed to restore gauge invariance in calculations where the 
numerator function of the 1-loop Feynman diagrams is computed in four dimensions. This 
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Figure 14: QCD Feynman rules (box on the top) and corresponding e-propagators and e-vertices. 
q represents the flow of the loop momentum. 
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possibility is particularly appealing in conjunction with schemes such as dimensional re- 
duction, where the use of particular classes of four-dimensional identities involving 75 is 
forbidden. 



It would also be interesting to generalize this approach to the Unitary gauge and 
beyond 1-loop. I leave these two issues to future investigations. 
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A. Prom the FDH scheme to the HV scheme 



In this appendix, I give the translation rules from the FDH scheme (or dimensional re- 
duction) reproduced by the Feynman rules given in sections ^ to |^ and the HV scheme 
of eq. (1.5). For QCD, once a 1-loop amplitude A^^^ has been computed in FDH, the 



corresponding HV result can be obtained with the help of the formula |]69| 



4(1) - /iW + 4(0) 9^ riV.. n„ 



^HV - ^FDH ■ -- ]^g^2 



+ -^>-f^ , (A.1) 



where A^^^ is the tree level result, Nc the number of colors, Uq the number of massless quarks 
and uq the number of massive quarks. For the electroweak model, once the renormalized 
1-loop amplitude has been determined in FDH, since the terms proportional to e in eq. ( |1.5| ) 
are separately gauge invariant, their total contribution can be completely reabsorbed by 
shifts of the renormalization constants. The HV result can then be obtained through the 
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replacements ^ 
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f,R 
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247r2s2 
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'l27r2s 



247r2 



1 



+ 



+ 



2 167r2 I 4s2c2 
e2 / 1 



^C2 ^2s2 



167r2 \ 4s2c2 



.2 Q) 



c2 c2 2s2 



(A.2) 



167r2 c2 



From eq. ( [A.2| ), the necessary shifts in the charge renormalization constant and in the sine 
and cosine of the weak mixing angle (relevant when using the on-shell scheme) can be 
determined from the equations 



1 / s 

5Ze = --{ 6Zaa + -SZzA 

2 \ c 



6c 
c 

6s 
s 



1 /6M^ 6Ml 



Ml 



c2 Sc 
s2 C 



(A.3) 



The rules in eq. (A.2) are easily derived from the explicit knowledge of the part proportional 



to Ahv in the 2-point functions listed in [47|. 



''I use the same notations and conventions of l6ql and assume a unit CKM matrix. 
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